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We show that a thermal reservoir can effectively act as a squeezed reservoir on atoms that are 
subject to energy-level modulation. For sufficiently last and strong modulation, for which the 
rotating-wave-approximation is broken, the resulting squeezing persists at long times. These effects 
are analyzed by a master equation that is valid beyond the rotating wave approximation. As an 
example we consider a two-level-atom in a cavity with Lorentzian linewidth, subject to sinusoidal 
energy modulation. A possible realization of these effects is discussed lor Rydberg atoms. 
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I. INTRODUCTION 

A two-level atom (TLA) that interacts with squeezed 
light fields (or squeezed reservoirs) may exhibit charac- 
teristic dynamics and fluorescence properties that depend 
on the strength and/or phase of the squeezing [TJ|2]. For 
example, different decay rates for the two atomic-dipole 
quadratures [3] and subnatural linewidth of resonance 
fluorescence peaks g], were predicted. Some of these 
features persist when the more realistic non-minimum- 
uncertainty state of thermal squeezed light is considered 
[5]. Recently it has also been shown that the error in- 
duced in coherent control of atoms by quantum fluctu- 
ations and light-atom entanglement may be either sup- 
pressed or enhanced when a squeezed light pulse is used 
for such control [6j. However, a challenging limitation on 
the observability of most of these phenomena is the rather 
poor spatial overlap between the dipole radiation pattern 
of the atom and available paraxial sources of squeezed 
light. For this reason, atomic interaction with squeezed 
light in a cavity has been considered [7J. 

A completely different approach towards observing 
such phenomena involves emulating the interaction of 
atoms with squeezed light, while the atom actually in- 
teracts with a reservoir in a thermal or vacuum state. 
This may be achieved by controlling accessible param- 
eters of the atomic system without affecting the state 
of the inaccessible reservoir, thereby creating an effec- 
tively squeezed reservoir. For example, in [5] two states 
of a laser-driven four-level-atom are shown to be coupled 
to an effectively squeezed reservoir whose properties are 
determined by the laser parameters. More proposals us- 
ing multilevel atoms can be found in [5H11). In [T^] it 
is shown that when a TLA is driven by a strong classi- 
cal field, while coupled to a thermal reservoir, it may be 
viewed as coupled to an effectively squeezed reservoir, as 
long as the Rabi frequency of the driving field is not too 
small compared to the spectral width of the reservoir. 

Here we show that modulations of the level-spacing of a 
TLA coupled to a thermal reservoir can lead to the same 
effect as its coupling to a squeezed reservoir. The key 
principle of our approach, that allows for effective squeez- 
ing, is the breakdown of the rotating-wave-approximation 
(RWA) of the field-atom interaction, resulting in an inter- 



action similar to that obtained for a squeezed reservoir. 
The RWA implies neglecting the counter-rotating terms 
in the system-reservoir coupling that oscillate at least as 
fast as the level-spacing frequency, w a , and is valid for 
times much larger than u)^ 1 [S]- However, at shorter 
times, the counter-rotating terms are responsible for the 
appearance of terms characteristic of a squeezed reservoir 
in the master equation of the system, although it is cou- 
pled to a thermal, non-squeezed, reservoir, as was shown 
in [14] for quantum Brownian motion (conversely, when 
adopting the RWA but considering the interaction with a 
squeezed field mode, the effect of squeezing is reminiscent 
of the existence of counter-rotating terms 0). The ques- 
tion is whether modulation of the level-spacing can help 
preserve the counter-rotating terms even at times much 
longer than cu^ 1 and thus lead to coupling to an effec- 
tively squeezed reservoir. Such modulations have been 
extensively studied for both short, non-Markovian, time 
scales and long, Markovian, ones, and shown to dras- 
tically modify atomic decay [15], decoherence [TB] and 
thermodynamics [17] . through either Zeno-like suppres- 
sion of the coupling to the reservoir or anti-Zeno-like en- 
hancement. 

We find that when the level-spacing modulations of 
a TLA that is coupled to a thermal reservoir, are suffi- 
ciently fast and strong to break the RWA, the practically 
achievable degree of squeezing does not exceed that of a 
"classically squeezed state" [I], namely, a state wherein 
the effective squeezing parameter of the reservoir, M, and 
its effective mode occupation at w a , N, satisfy \M\ < N. 
In such a state the quadrature noise of the reservoir can- 
not be less than that of the vacuum [TJ [TS] . Emulating 
quantum squeezing, so that N < \M\ < y/N(N + 1) [1], 
is possible, in principle, by this method, but rather im- 
practical using present-day technology. In order to pre- 
serve the counter-rotating terms at long times, i.e. make 
them non-oscillatory in time, the required modulation 
strength and frequency is found to be of the order of 
uj a , which makes this scheme suitable for the microwave 
regime rather than the optical one. 

The paper is organized as follows. In Sec. II we 
provide simple arguments that explain the relation be- 
tween counter-rotating terms and effective squeezing, and 
the effect of modulations. Sec. Ill is devoted to the 
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derivation of a generalized non-Markovian master equa- 
tion (ME) for a TLA in a thermal reservoir, in the pres- 
ence of level-spacing modulations and without invoking 
the RWA. It contains our main general result, presented 
m Eqs. ( |19|20P , i.e. the emergence of effective squeezing 
terms in the ME for the system density operator. The 
general analysis is illustrated in Sec. IV, for the case 
of an atom coupled to a cavity with a Lorentzian spec- 
trum. In Sec. V we discuss the analytically solvable ex- 
ample of sinusoidal modulation and explicitly derive the 
corresponding squeezing parameter. Sec. VI includes a 
discussion of possible observable effects, such as atomic- 
dipole dephasing and fluorescence, in the context of our 
modulation scheme. Finally, Sec. VII presents relevant 
experimental considerations and an example of a possi- 
ble Rydberg-atom realization of such effects. The con- 
clusions are presented in Sec. VIII. 



II. WHY COUNTER-ROTATING TERMS 
INDUCE EFFECTIVE SQUEEZING 

In this section we wish to explain, in simple terms, our 
main idea for engineering an effectively squeezed ther- 
mal reservoir. We begin by considering the interaction 
between a TLA with energy huj a and a broadband elec- 
tromagnetic field with resonant carrier frequency uj ai 



H = hd (<7+ + <t_) (Ee~ iUat + E ] c 



(1) 



where E is the positive-frequency envelope of the field, d 
is the dipole matrix element and a± are the TLA opera- 
tors. In the interaction picture we obtain 

Hi = M (a+E + a -if 1 + a+E^e 12 ^* + a-Ee~ i2 " at 

(2) 

where E(t) is time-dependent. The last two terms are 
usually neglected for times t>%' in the RWA. Without 
invoking the RWA, this Hamiltonian can be written in the 
form, 



(3) 



F(t) = d\E(t) + £t(i)e i2w «' 



where F(t) may be recognized as a Langevin-like force, 
due to the reservoir, acting on the atomic-dipole. In the 
Markov approximation, the Langevin force is assumed 
delta-correlated in time, hence the normal and anomalous 
correlations of the force then determine the parameters 
M and N, respectively, as follows: 



(F\t)F(t + T)) = 1 NS(r), 
(F(t)F(t + r)) = jM5(t), 



(4) 



where 7 is the decay rate of the TLA population due to 
the reservoir. Let us now compare the origin of M and 



N in the cases of (a) a real squeezed reservoir with RWA, 
(b) thermal reservoir without RWA and (c) the case of 
modulations. 

Case a: squeezed reservoir + RWA. — assume that the 
broadband field E(t) has squeezed field correlations, i.e. 



(E(t)EHt 
(E\t)E(t 
(E(t)E(t 



■r)) 
■r)) 
r)) 



(n+l)S(r), 

nS(r), 

m5{j). 



(5) 



By taking the RWA in Eq. (J3J) we recognize that F — dE, 
and comparing Eq. Q to Eq. ^ we obtain M = d 2 m/j, 
N = d 2 n/-f, i.e. the effective reservoir, F, is squeezed 
(M ^ 0) since the "real" reservoir, E, is squeezed (m 7^ 
0). 

Case b: thermal reservoir without RWA. — consider 
now that the broadband field E{t) is thermal, namely it 
has only normal correlations, n^O and m — in Eq. ([5|. 
However, if we assume that the relevant timescales we are 
interested in satisfy t lj~ , then we may set e l2uJat « 1 
in Eq. ^ and the effective Langevin force becomes F = 
d(E + E^). Plugging this force into Eq. (|4| and using 
Eq. ([5} with m = 0, yields M = N = d?(2n + 
This shows that (1) the system is effectively coupled to 
a squeezed reservoir due to the existence of the counter- 
rotating terms, and (2) the achievable squeezing by solely 
breaking the RWA is that of a " classical squeezed state" , 
i.e. the magnitude of the anomalous correlation M does 
not exceed the normal correlation N . 

Case c: thermal reservoir + modulation, for long 
times. — here we take again a thermal reservoir, n 7^ 
and m = 0, and take the long-time limit t ^> uj~ 1 only 
when modulations are considered. The modulated TLA 
Hamiltonian is taken to be 



H A (t) 



1 



: ha. 



S(t)} 



(G) 



Then, in the interaction picture &-j(t) = <r_e lUat e*(t), 
with 



(7) 



where we assumed that e(t) can be written as a Fourier 
series, and let us take its period to be 27r/w a , i.e. v q — 
qui a with integer q. The Langevin force, Eq. ([3]), be- 
comes, 



F(t) 



E{t)e 



iqoj a t 



£t( t y(2+<r)toa 



(8) 



From the above equation it is evident that in order to ob- 
tain squeezing similar to that in case b, it is sufficient that 
the terms q — and q = — 2 exist. Then, at long times 
t w" 1 , at which the oscillatory terms e t2uJat vanish, 
we obtain F(t) w d(e E + e_ 2 E), with M = e e_ 2 (2n + 
l){d 2 h) and N = [|e | 2 n + |e_ 2 | 2 (n + I)](d 2 / 7 ). For 
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|eo| = | £-2 1 we recover case b where classical squeez- 
ing is achievable. Quantum squeezing however, namely 
\M\ > N, is practically more difficult to obtain, though 
possible in principle, as explained in Appendix A. At this 
point it also becomes clearer why strong and fast modula- 
tion 5(t) is required: the existence of squeezing is enabled 
by the q = —2 term with frequency — 2w a . Then J dt'S(t') 
has to contain a frequency of order u> a , e.g. sin(w a i), such 
that the derivative 8{t) includes a term with magnitude 
and frequency w a , e.g. w a cos(w a i). This summarizes our 
basic idea of how to create an effective squeezed reservoir 
by system energy-modulation. 



III. GENERALIZED MASTER EQUATION 

To provide a more rigorous account of such effects, 
we derive a master equation (ME) for the TLA density 
operator, in the presence of electromagnetic field modes 
in a thermal state, when the TLA level-spacing undergoes 
modulations. The system+reservoir Hamiltonian reads, 



term in (12), and, more specifically, on three representa- 



H(t) = H A (t)+H F 



Haf = 



E 



hgkip- + cr+)(<4 + a k ), 



(9) 



where k are indices of different field modes with corre- 
sponding destruction operators a k , frequencies u>k and 
dipole couplings g k , and Ha is the one in Eq. In the 
interaction picture with respect to Ha + Hp, the Hamil- 
tonian becomes 



Hj(f) = h\a-Fl{t)e*{t)+d_F a e*{t)+Yi.c, 



Fa{t) 



fe 
fe 



"a) 



(10) 



with e(t) from Eq. §f§. Note that the tilded operators 
F a (t) , (t) are the ones that are neglected in the RWA 
without modulation. In the Born approximation, namely, 
when neglecting the system-reservoir correlations, the 
ME for the TLA density operator p is [19] , 



m = 



dt'tTp([H I (t),[H I (t'),p(t')p F }}), (11) 



where pp is the stationary density operator of the reser- 
voir (the field). The double commutator inside the trace 
over field degrees of freedom contains four terms, 

H I (t)H I (t')p(t')p F - Hj^pit^ppHjit') 
-H^p^ppH^t) + p{t')p F H I (t')H I {t), (12) 

each of which contains 16 terms. In order to illustrate the 
main points of the derivation, let us focus on the second 



tive terms out of its 16 terms 

a + p(t')a_F a (t)p F Ft(t')e(t)s*(t'), 
& + p(t')a+F a (t)ppF a (t')s(t)e(t'), 
& + p(t')a + F^t)p F F a (t')s(t)e(t'). (13) 

The first two terms do not contain tilded operators, so 
they exist also in the RWA without modulations. When 
the trace over the field is taken, it can be seen that the 
first term includes the normal correlator (F^(t')F a (t)) = 
tr p(F^(t')F a (t)pF), which exists for a thermal reservoir. 
This term contributes to the a + p&- term in the ME, 
typical of a thermal reservoir jTSJ HH] . The second term 
involves the anomalous correlator (F a (t')F a (t)) which ex- 
ists only for a squeezed reservoir and contributes to the 
squeezing term cr + pa + in the ME |18j . However, in our 
case we assume a thermal, non-squeezed, reservoir, and 
this correlator vanishes, so that indeed, for RWA terms, 
no squeezing is expected when a thermal reservoir is con- 
sidered. The third term in ( 13 ) is the most interesting 



for our purposes. It does not exist in the RWA as it con- 
tains a tilded operator, yet, if the modulations preserve 
it at long times, it contributes to the squeezing terms of 
the ME, a + pa + . We calculate its associated correlator 
(F a (t')F^(t)) using Eq. Iiol, 



(^(^l(i)) = EE«t" 



-i(w h -u a )t'i(ui' k -t 



k k' 



42 ^* / du J D{uj)\g(uj)\ 2 e- t ^ +UJ ^ it - t " ) [n(Lj) + 1]. (14) 



Here we used {d\dk>) = Skk'n>k for a thermal state of the 
field, where n k is the Planck distribution for frequency 
Wfe, and we defined the density of field modes by ^ fe — > 
J dujD(uj). We further define the coupling spectrum of 
the reservoir, 



GoM = D(u)\g{u)t 



(15) 



Reexamining the expressions in (13) and the ME flllk, 



we need to multiply this correlator by the modulation 
functions and integrate over time, 



£ a >e 



duG (u)[l + n(u)] / dt'e- i ^ +u " +v ^ t - t '\ 



(16) 



where the definition of e(t) from Q was used. Here p(t') 
was taken out of the integral by assuming that t is much 
shorter than the typical timescale for changes in the sys- 
tem, i.e. p{t') ~ p(t). In cases considered later, when the 
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Markov approximation is taken and the resulting ME has 
time-independent coefficients, we may view this short t 
as a coarse-graining time and the validity of the ME is 
then extended to long times |13j . In Appendix B we show 
that the term in ( 16 ) can be written as 



Y,Y. £ i £ <i' el(Uq+Vq ' +2Ua)t 

q q> 

"1 



du)5t{ui) 



7„(w -Ul a - Vq) - lA„(w -U) a - Vq) 



(17) 



where S t (uj) = ^ sm ^" t - ) is a sine function that approaches 
a delta function as t gets larger, and 



7„(cj) = 2ttGoH[1 + n(w)] 
A„(cj) = P J 



<L/ G V)[l + n( t S)] t (18) 



uj' — UJ 



with P denoting the principal value. Repeating the pro- 
cedure shown in Eqs. (14|16|17l for all terms of Eqs. 



( 11|12 ), we find the following generalized non-Mar kovian 
ME, 

1 7 

P = -ih A& z>P\ + o N[2a + pd- - po-& + - a-a+p] 



-{N + l)[2<r_p(7 + - pfr+a- - a+a-p] 
r )Mo + p<j + — jM*a—pa—. 



(19) 



The first term of the ME is a Hamiltonian term contain- 
ing the TLA energy correction HA, and would not be of 
interest in the following. The second and third terms 
describe absorption and emission from and to the pho- 
ton field, respectively, typical of damping by a thermal 
reservoir, where 7 is the decay rate and N the reservoir 
population. 

The last two terms, which are the most interesting for 
our discussion, appear to be squeezed-reservoir damp- 
ing terms. Although we assumed a thermal reservoir, 
they exist due to the modulation. The ME coefficients, 
A, j,N are real (whereas M is complex) and are gener- 
ally time-dependent and affected by modulation: 



A 



-yN 



j(N+l) 



7M 



'/</' 

-i 



qq' 



dujS t (uj) [Ar(u + w n + v q ) - A t (uj — w a — v q )] 



duS t (u}) - [7 T (w + uj a + v q ) + 7 T (w - w a - v q )] , 

1 « 



r i Z* 00 



+z [sqE*^-^ - e^e-^-^] / dw«t(a;)A r (w + oj a + v q ), 



— OO 
OO 



duSt (u>) - 7t (w + w a + ^ 9 ) 



^( w )[-7 T (a; + a; +^) + -7 T (. 

w 

+iA T (uj + L> a + Vq) - iA T (uj -w a - v q )}, 



LO -UJ a ~ V q ) 



(20) 



where we used the temperature-dependent response of long times, when (5 t (w) is treated as Dirac delta, the mod- 
ulation enables a scan through the frequency response of 
the reservoir 7r(w) + iAt(w) and thus change the effec- 



tive reservoir 

7t(w) = 2ttGt{u) 



A T (w) = P duj 



,G T (u/) 



tive coupling to it. This was also concluded in |16j . where 
modulations were shown to be useful for decoherence con- 
trol. As for the squeezing, 7M, let us first consider it 

GtM = G {u)[l + n(u)] + G {-u)n(-u). (21) without modulation, i.e. when e = 1, v Q = 0, = 0. 

Then 7M is fast-oscillating and negligible for t ^> u a 1 , 



The effect of the modulation is clear from Eq. ( 20 ) . At 
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as it should in the RWA. However, by choosing the right 
modulation frequencies v q this term may become impor- 
tant at long times as well, as was explained in Sec. II 
above. 

The above ME may be viewed as a generalization of 
previous treatments [El [16] : (1) here the principal- value 
terms A-r(u; a + v q ) are not neglected, and can in fact 
become important for reservoir spectra Gt{oj) which are 
asymmetric around u a + v q . (2) The ME is written in 
operator from, rather than as a Bloch-equation for the 
TLA density matrix elements, which makes it easy to 
generalize to other systems, e.g. an harmonic oscillator 
(see Sec. VIII). 



IV. LORENTZIAN RESERVOIR: ATOM IN A 
CAVITY 



y T /Yc, A T /Yc 

1000 




FIG. 1: (color online). Temperature dependent response of 
the cavity reservo ir, 7t(<^) (blue solid line) and At(w) (red 
dashed line) , Eq. ( 25 I . Here j c is the atomic decay rate to the 
cavity reservoir in the absence of modulation, and u> a = 10k 
is taken. 



Up to this point we have not specified what are the 
field modes which make up the reservoir. Let us now 
consider an atom in a resonant cavity in the " bad-cavity" 
limit. The interaction of the atom with a cavity-mode 
with envelope a can be written in the interaction picture 
as 



Ei = ng[a(t)e- iUat + a f (^^(a+e 1 ^ + 



— iuj a t\ 
)■> 

(22) 

with dipolar coupling g, so that the Langevin force acting 
on the atom in analogy with Eq. (10) is F a — ga{t). 



The cavity-mode is damped by the coupling to outside 
modes and its correlation functions decay exponentially 
(see Appendix C for details), 



(a{t)a){t + T)) 
(at (t)a(t + t)) 



l)e 



-fM 



n a e-fl T l, 



(23) 



where n a — n(uj a ) and k is the width of the cavity-mode 
due to the coupling to outside modes, assumed here to 
satisfy k ^> g (bad-cavity limit) and k <C ui a (see Ap- 
pendix C). 



A. Lorentzian response 

From the Fourier transform of the correlation func- 
tions above we can deduce the Lorentzian spectrum of 
the cavity-mode, which makes up the reservoir for the 
atom, 



G (w)[n(w) + 1] -> g 2 {n a + l) — 



( f)2 + (w _ Wa) 2 



(f)2 + ( w _ WB )S 



(24) 



Using this spectrum in Eq. (21 1 we find the temperature- 
ionse of 1 

K + i) 



dependent response of the reservoir 

1 



7t(w) 
A T (w) 



7c 



7c 



(«o + 1) 



l+4(^^) 2 

K 2 (uJ a - LP) 
4(u! a — Uj) 2 K + K 2 



1 



K 2 (w a + Uj) 
a 4(w a + U>) 2 K + K 2 



7 C = 27rG (w a ) = 2ir 



4 5 2 



(25) 



where j c is the TLA damping rate to the cavity reservoir 
without modulation. The real and imaginary response 
functions obey the Kramers-Kronig relation and are plot- 
ted in Fig. 1. 



B. Long-time approximation 



Going back to Eq. (20) for the ME coefficients, we 



recall that the sine function 5t{ui) is of width 1/t. The 
typical scale of variations of the response functions ( 25 ) 
is K, so by assuming t ^> k _1 we may take <5t(w) f=s SJuj) 
in the integrals over the response functions. This is the 
well-known Markov approximation, which means that we 
coarse-grain over timescales of the order k" 1 , such that 
all the observed phenomena should be slower. Moreover, 
we are also interested in timescales much longer than 
a;" 1 , as in the RWA. In Sec. II we saw that in order to 
obtain effective reservoir squeezing for t ^> uj^ 1 , the mod- 
ulation frequencies v q should be at least of the order u a . 
This means that we should keep only the non-oscillatory 
terms in the ME coefficients, i.e. v q — v q > in 7 AT, j(N+l) 
and v q + v q > + 2u a = in jM. To sum up, taking the 



approximation t 3> 



1 , we obtain for these ME 
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coefficients, 

7 iV 

7(JV + 1) 
7M 



E \ £ q?lT{-U a ~ Vq), 
1 

J2\£q\ 2 lT(Ua + Vq), 



qq' 



+2bJ a )t 



[)p T {Ua + Vq) + ]p T {~U a - Vq) 

+iA T (uj a + v q ) - iA T (-uj a - v q ) 



(26) 



where J2 qq ' denotes the sum with the constraint v q 



20J n 



0. Here the above coefficients are time- 



independent and the ME, Eq. ( 19 ), becomes a Markovian 



ME of a system damped by an effective thermal squeezed 
reservoir, which is valid at long times. 



V. SINUSOIDAL MODULATION 

In order to illustrate our scheme for engineering an 
effective squeezed reservoir for the TLA, let us take the 
analytically solvable example of a sinusoidal modulation 
of the TLA level spacing, 



5(t) = zmuj a [l — aisx(rnu} a t)], 



(27) 



where z and m are positive parameters. Introducing the 
above modulation in Eq. ([7]), and using the identity 
e izcos0 _ Y^ q x L_ 0o %i J 9 (z)e l9 ^with integer q, we obtain 

e- iz i«J q (z), 



v q = (z + q)muj a . 



(28) 



The constraint v q + v q i + 2us a = imposed on the 7M 
term in the long-time approximation now becomes q' — 
— (2/m) — 2z — q, and since q', q are integers, we get the 
following constraints on the parameters z and m, 



integer 



integer 



(29) 



The expressions for the coefficients from Eq. ( 26 ) now 
become, 



OO 

E 1 

q—~oo 



Jq( z hT(-u a [l + mz + mq]), 



7(JV+1) 



Jqizhriuall + mz + mq}), 



OO 

E 

q— — 00 



e -«'(-l)s+»J,(z)J. 



-q—2z- 



.(z)x 



[I 7T ( Wa [l + mz + m(? ]) + I 7T (_ Wa [i _|_ mz _|_ jjiq} ) 

+iAr(a; a [l + mz + mq]) — iAx(— cj a [l + mz + mq])]. 

(30) 
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FIG. 2: (color online). Effective reservoir squeezing 7IMJ as 
a function of modulation parameters z and m, Eqs. | 27|29| |. 
Taking oj a = 10k and n a = 10 in Eqs. | |25|30[ ), the difference 
between and M is plotted in (a), (b) and (c) and apparently 
cannot go lower than the "classical squeezing" bound \M\ < 
N, as argued in Sec. II. (a) m=2, (b) m=l, (c) z=l. (d) The 
ratio \M\/(N + 0.5), that determines the difference between 
the two atomic-dipole quadrature dephasing rates, is plotted 
for m = 2. 



Considering the Lorentzian reservoir from ( 25 ) with uj a = 
10k and n a = 10 3 , let us choose for example z = 1 and 
m = 2. It is enough to take q from —20 to 20 in order for 
the sums in Eq. (30) to converge, and we obtain, 7^ = 
207.49 7c , 7(JV + 1) = 207.67 7c , 7 M = (92.08 - i42.14) 7c 
and | 7 M| = 101.267 c . This demonstrates an appreciable 
effect due to modulations, as \M\ becomes about half of 
N. On the other hand, we also verified that for z = we 
return to the RWA results without modulation, namely, 
lN = 7 c n a , 7 (A + 1) = 7c (n a + 1) and 7M = 0. In 
Fig. 2 we plot, with the same ui a , n a and as a func- 
tion of the parameters z and m, the difference between 
"f\M\ and 7 7V. This quantity determines the increase 
of the atomic-dipole dephasing rate and broadening of 
the fluorescence line shape, as discussed in Sec. VI be- 
low. As this difference becomes smaller, the increase 
and broadening become smaller compared to the non- 
squeezed reservoir case. Figs. 2a and 2b present this 
quantity as a function of z for fixed m values, m = 2 
and m = 1 respectively, whereas in Fig. 2c it is plotted 
as a function of m for z = 1. It appears that m = 2 is 
optimal, and in the range of z values considered, z = 9.5 
provides the smallest result, 7A — j\M\ — 0.128 7c with 
7 iV = 0.683 7c . This suggests that indeed \M\ does not 
exceed N, i.e. we get "classical squeezing" by breaking 
the RWA, as argued in Sec. II. 



VI. POSSIBLE OBSERVABLE EFFECTS 

In this section we wish to describe possible observable 
manifestations of the effect of modulations and effective 
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squeezing. These include two distinct dephasing rates 
of the atomic-dipole quadratures, fluorescence and reso- 
nance fluorescence. 



A. Atomic dipole dephasing 



From the master equation (ME), Eq. (19 1, we ob- 
tain the equations of motion for expectation values of 
the atomic operators, 



(&+) = -7(^+0 (a + )- 7 M*(a. 
(& z ) = - 7 (2iV + l)(^)- 7 , 



(31) 



and (&_) = (&+)*. Denoting M = \M\e i2ip , we use the 
squeezing phase tp to define the atomic-dipole quadra- 
tures as a x — e ltp (T^ 
ie~ ltp (j-, and obtain 



e lv a- and a y = 



(*x) = -7 [N - \M\ + -J (a x ) = - lx (a x ), 

(& v ) = -7 U+ \M\ + (a y ) ee - 7y (a y ). (32) 

The above result, shown previously by Gardiner [3], 
reveals that the two different quadratures may decay 
with very different rates when the reservoir is squeezed, 
\M\ 7^ 0. In our case, \M\ does not exceed N, so that the 
slower rate 7^ is bounded by the vacuum rate 7/2. Note 
that the effect of modulation on the quadratures' dynam- 
ics is twofold: not only does it induce squeezing and hence 
two distinct decay rates, but it may also change N and 7 
from their thermal-reservoir values without modulation, 
n a and 7 C respectively. In cases where the dynamics can 
be measured, they may reveal both effects in the most di- 



t|m|| 



for m 



rect way. Fig. 2d portraits the ratio ^ N+0 5) 
and as a function of z. As can be seen in Eq. (32 1, this 



ratio determines the difference between the two decay 
rates. For instance, when 2 = 1 we obtain "f x = 1067 c , 
7 a = 3097c, whereas for z = 3 we get j x = 32, 7 y = 195, 
and the corresponding ratios for these z — 1,3 cases are 



7|M|| 
7(W+0.5) 



0.49,0.72. 



B. Atom fluorescence 



Making use of Eq. ( 32 1 and the quantum regres- 



sion theorem [T5], it is easy to obtain the correlation 
function C(t) and the corresponding spectrum S(uj) = 
- dte~ luJt C(t) for the atom in steady state, 



C(t) = (&+(0)*-{t)) 



1 TV 

2 2N + 1 



.-7x1*1 



-7„l*l 



1 



S(u) = - 



N 



tt2N + 1 



ly 



% + w 



(33) 
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FIG. 3: (color online) . Atom fluorescence spectrum, Eq. (33 1 



for modulation parameters z = 3, m = 2, and with cu a = 10k, 
n a — 10 3 . (a) the spectrum under the influence of modulation 
(blue solid line) is compared to that in the absence of mod- 
ulation (red dashed line, multiplied here by factor 5). The 
modulated case has a much narrower spectrum, (b) the spec- 
trum in the modulated case (blue solid line) compared with a 
Lorentzian of width 7(^ + 0.5) (red dashed line). It is appar- 
ent that the shape of the spectrum in the modulated case is 
not that of a single Lorentzian, as indeed can be seen in Eq. 
d33]. 



where ui here is shifted by uj a , i.e. the spectrum is actually 
centered around uj a . This spectrum is proportional to 
that of the fluorescent light emitted by the atom and 
it reveals both effects of the modulation. First of all, 
even if squeezing was absent and the spectrum would 
still be a Lorentzian of width ^(N + 0.5), the modulation 
changes it from the unmodulated thermal value, (n a + 
0.5)7 C = 1000.57 c . Perhaps the more interesting effect 
however, is that of squeezing 3J, namely the spectrum is 
actually not a single Lorentzian, rather it is a sum of two 
Lorentzians with different widths 7^ and jy. As in the 
dipole dephasing case, the effect becomes more apparent 
as the two widths become more distinct. In Fig. 3 we 
present, for the case m — 2, z — 3, the comparison of the 
fluorescence spectrum with and without modulation. We 
also show the non-Lorentzian shape of the spectrum by 
comparing it to a Lorentzian of width 7(iV + 0.5) in the 
same case. 



C. Resonance fluorescence 

The resonance fluorescence spectrum of an atom in a 
squeezed reservoir was treated in [4] and extended to the 
case of a thermal squeezed reservoir in [5J. The atom 
is driven by a resonant strong (classical) field with Rabi 
frequency fie 10 , and the resulting fluorescence spectrum 
is given by [U [5] , 



S [a.u.] 



S(oj) oc 



70 



70 



(w - w a ) 2 + 7</> 2 

r,* 




(w - w a + n) 2 + r 2 ( w - w - ft) 2 + ry 
i 



7 [N + |M|cos<£ + 



-1500-1000 -500 

(b) S [a.u.] 

0.03 



500 iooo i5o(r yc 



1 



r = (n-^\m\ coa<r , }) . 



(34) 



with <p — 2(8 — ip). The spectrum is comprised of the 
three well known Lorentzian peaks centered at lo ci and 
cj a ± 51 |20| . but here it is modified by squeezing. In fact, 
we obtain a phase-dependent phenomenon, namely the 
widths of the fluorescence peaks, 7^ and T^, change when 
the phase of the strong field is varied. The phase depen- 
dence is only due to squeezing and provides a very good 
evidence of the effective squeezing induced by the modu- 
lation. In Fig. 4 we plot the central peak of the spectrum, 
a Lorentzian of width 7^, for z = 1, 3 and different values 
of <f). Taking the case z = 1 for example, and beginning 
with <f> = 7r/2, \M\ cos cf> — and the only effect of the 
modulation is the modification of the reservoir modes 
occupancy N with respect to the thermal one n a . This 
results in a narrower peak, 7^ = -f(N + 0.5) = 2087 c in- 
stead of 1000.57 c without modulation. When <f> is varied 
to 0, 70 gets a positive contribution from the squeezing 
term \M\ cos</f> = \M\ and as a result the peak broadens 
to 70 = 7^ = 3097 c . The narrowest peak is obtained for 
(j> = 7r for which 7^ = j x — 1067 c . 



VII. EXPERIMENTAL CONSIDERATIONS 

We need to find a dipole transition, whose resonant en- 
ergy uj a can be controlled such that the two levels are still 
well separated from others. Moreover, in the sinusoidal 
modulation considered in Sec. V, the amplitude and fre- 
quency must exceed u> a . This leads us to consider two 
almost degenerate states whose energy splitting may be 
controlled by oscillating fields in the microwave regime. 
Since the effects discussed here involve coupling to radi- 
ation, it is also preferable that the dipole moment of the 
transition be large, as for Rydberg levels. 

Let us now discuss the timescale of the above phe- 
nomena. Both the dipole dynamics and the fluorescence 



spectral widths are appreciable for t > y x 



1 hi 



The 



two rates ^ x y become distinct when there is squeezing, 
and in the cases discussed here, are of the order of O.I7T 
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FIG. 4: (color online). Resonance fluorescence spectrum un- 
der the influence of modulation and effective squeezing, as a 
function of the relative driving field phase (f>, for <f> = tt (blue 
solid line), (f> — n/2 (red dashed line) and <f> = (black dotted 
line). The modulation parameters are: (a) z — 1, m = 2, and 
(b) z — 3, m — 2. As in Figs. 2 and 3, uj a = Wk and n a = 10 3 
are taken. 



to 7r, with 7^ = 7 c w a rj j c T/(Hcj a ), where T is the 
photon-reservoir temperature in energy units, and is 
the corresponding width without modulation. Hence, we 
wish to observe phenomena at times t > 7^ 1 - 

Re call ing the atomic damping in the cavity, 7 C from 

(25), we have j T = 



Eq 



T 



with q = d 



2e hV ' 

Here d is the dipole matrix element of the TLA transition 
and V ~ is the cavity mode volume, taken here to be 
of the order of the transition wavelength cubed. We thus 
obtain 



It 



3 T 
2tt ^° Hk ' 



(35) 



where 70 



is the atomic transition spontaneous 
The above expression shows 



37T€q?IC 3 

emission rate in free space, 
that the considered effects grow with the temperature 
divided by the cavity linewidth, so that high temper- 
atures and narrow-linewidth cavities (still in the bad- 
cavity limit) are preferable. 

Finally, let us recall the conditions of validity for some 
of the approximations taken in our analysis. The deriva- 
tion of the reservoir spectrum in the bad-cavity limit re- 
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quires u> a 3> k and k S> g, 70 • The long-time limit in Sec. 
IV B is valid for ( > k" 1 , i.e. in our case we demand 

Example: Rydberg levels. — Consider the circular Ry- 
dberg level \e) = \n = 51,1 = 50, m = 50). By dipole 
selection rules, it can only decay to the lower circular 
level \s) = \n = 50,1 = 49, m = 49) or to the almost 
degenerate level \g) = \n = 51,1 = 49, m = 49). We 
choose |e) and \g) as our excited and ground states, re- 
spectively, and by a static magnetic field can induce an 
energy (frequency) shift, of say, 1GHz (the s O e transi- 
tion frequency is 51.1GHz). Another magnetic field, this 
time RF, is used for the modulation. It oscillates at a 
frequency of 2GHz and is strong enough to induce an en- 
ergy shift of a few GHz. For example, a scheme for atom 
fluorescence may consist of first preparing the atom in 
state |e), then applying modulations and measuring the 
fluorescence spectrum. Since the atom may also decay 
to the lower circular state \s), the spectrum may have an 
additional peak at 51GHz, but the interesting peak for 
our purposes is the one at oj a /(27r) = 1GHz. 

To evaluate the typical timescale of such an exper- 
iment, we calculated the dipole matrix element of the 
above g •<-> e transition by assuming that for such high 
n numbers the electronic wavefunctions are hydrogenic 
(see Appendix D). This is such, since in these large n's 
the last electron is far from the rest of the atom and ef- 
fectively "sees" a unit positive charge like the Hydrogen 
nucleus. We obtained, e.g. for the x component of the 
dipole, d = 382ea , where e is the electron charge and 
ao the Bohr radius. The spontaneous emission rate of 
the transition then becomes 70 ~ 1.1 x 10 _5 Hz, and by 
taking k = 2x 10 Hz, we obtain, at room temperature 
T = 300K, 7 T w 2162Hz. This means that the exper- 
iment should last 0.4ms to 4ms. We also verified that 
the conditions of validity mentioned above are satisfied. 
The effect can be further enhanced by effectively increas- 
ing the reservoir temperature by shining on the atom hot 
thermal radiation, or the radiation from any other strong 
incoherent source, which is wideband around u> a . 



VIII. CONCLUSIONS 

This study used the apparent equivalence between 
counter-rotating terms in the reservoir-atom interaction 
and the relaxation of an atom coupled to a squeezed 
reservoir, to show how to engineer a thermal-reservoir 
so that it emulates a squeezed-reservoir. The key point 
that allows us to achieve this goal is the breaking of the 
rotating-wave approximation even at long times, by fast 
and strong system energy-modulation. To this end, we 
derived a general master equation, without invoking the 
Markov and rotating-wave approximations, so as to ren- 
der the effect of modulation in clear form. The master 
equation assumes a two- level system, yet, it is easily gen- 
eralized to any system-reservoir interaction of the form 
XF, where A" is a system operator with matrix elements 



between different energy states of the system and F is a 
reservoir operator with a continuous spectrum: e.g., for 
a system of a harmonic oscillator with lowering opera- 
tor a, X = a + fiT, and in the master equation cr1,<7+ 
are replaced by a, (with a modification for the en- 
ergy correction term A) . Extension of the treatment to 
energy-level modulation in multilevel or multipartite sys- 
tems described by large angular momenta |22j is expected 
to yield qualitatively new features. 

In order to illustrate our scheme for squeezed reservoir 
engineering, we considered the case of an atom inside 
a cavity, namely a Lorentzian reservoir, in the Markov 
regime, and applied sinusoidal modulations. We dis- 
cussed possible measurable effects in atomic-dipole de- 
phasing and fluorescence, along with experimental con- 
siderations and a possible realization of a system com- 
prised of two atomic Rydberg levels. 
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Appendix A: Emulating a quantum squeezed 
reservoir 

In Sec. II, we found the resulting effective reservoir 
parameters M and N for case c, where modulations are 
considered and only the £0 and e_2 terms exist. Looking 
at their difference, 

N - \M\ cx (|e | - |e_ 2 |) [n(|e | - IM) - IM] - (Al) 

we can see that quantum squeezing, i.e. \M\ that satisfies 
N < \M\ < y/N(N + 1), becomes largest for n = 0, i.e. 
for zero temperature, when the original reservoir has no 
thermal excitations. This is indeed what we observe for 
the Lorentzian reservoir in Sec. IV and the modulation 
from Sec. V. Taking n a — 0, co a — 10k and modulation 
parameters m = 2 and z — 1, nearly perfect quantum 
squeezing of M = 0.983 v/A(A + 1) is achieved, with 
j x = 0.0537 c an d ly = 0. 1547 c . The problem remains to 
find a realization using present-day technology. As men- 
tioned before, since fast and strong modulation is more 
suitable for microwave frequencies, and since here we saw 
that sufficiently low temperatures are needed for n„ « 1, 
it is natural to consider superconducting qubits in circuit 
QED cavities [53]. Then for w a = 2x 10 5 k = 2?r x 2GHz 
and typical dipole matrix elements d = 10 4 eao we find 
7c ~ 5.7KHz. The duration of the experiment is set by 
7^7 1 3.28ms, which is incompatible with the super- 
conducting qubit dephasing time which currently ranges 
between 1/j.s and 10/is. Future technology may allow em- 
ulating a quantum squeezed reservoir by the method de- 
scribed in this paper. However at present this appears to 
be very challenging. 
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Appendix B: Note on the derivation of the master 
equation 



Here we would like to show how to get equation (17) 



from Eq. (16 1. Beginning with the time integration over 
if and denoting v = —uj a — v q , we change variables to 
t = t — t' and get 



dt , e -i{ u -v){t-t') = / dTG(T)e~ i{u -" )T , (Bl) 



with Q(x) the Heaviside step function. By contour inte- 
gration one can obtain the relation, 



6(r)e' 



-i(o?— u)r 



lim 



1 



doj'- 



r;->-o+ 2iri J_ 00 uj' + if] — (uj — v) 
| ' (B2) 

Inserting Eqs. ( |B1|B2 ) into the double integral in Eq. 
( 16 1 (denoted here /), we get 



I 



rt e -iuj T 

duJ / dr—r lim / duj- 



-G n (u) 



_ t 2iri 7J-S.0+ J_ 00 uj' + ir] — (uj — v)' 

(B3) 

with G n ((jj) = Go(uj)[ti(uj) + 1]. Using the relation 

lim ~* r = mS{u-(u/+v))+P 1 ; 

t?— *-o+ uj + irj — (uj — v) uj — (uj' + v) 

(B4) 

where P denotes the principal value under integration, 
and defining the sine function 27r<5 t (u/) = J*^ t dTe~ lu r , 
we finally obtain, 



duj'S t (uj') 



iP I ,L: Gn ^ ; - irG n (uj' + v) 



uj + (uj' — v) 



(B5) 



Noting Eq. (18 1, this is just the expression in Eq. (17) 



with uj replaced by uj' . 



Appendix C: The spectrum of a cavity reservoir 



Our aim here is to derive the correlation functions in 



Eq. (23), for a cavity-mode in the bad-cavity limit. We 



begin with the model Hamiltonian H = Hs + Hr + H$r 
where the system (Hs) is the atom and cavity mode (a), 
the reservoir (Hr) is the electromagnetic modes that the 
cavity is leaking to (bj) and the perpendicular modes that 
the atom decay to (r;), and their interaction is described 
by H S r, 



Hr 



huj a a) a + 0.5%a z uj a + ^g(d + a^)(<r+ 



+&-), 



3 I 

it HA. —i— hh 



H S r = hJ2vj(a + d 1i )(b ] +St) 

3 

+h*^2gi(°- + v+) (n + 



(Cl) 



By writing the Heisenberg equations for d and bj and 
inserting the later into the former we obtain an equation 
for the envelope of the cavity mode a(t) — d(t)e lUat , 

a = - V |^| 2 f dt'a^e-^-^-^ + F(t) - igft, 

3 J ° 

F(t) = -iJ2vA(0)e- i ^-^ t , (C2) 



where fi = & + + cr_ and F is the Langevin force induced 
on the cavity mode by the outside modes. As in stan- 
dard Heisenberg-Langevin theory [3T], we now take the 
Markov approximation, i.e. we assume that the typi- 
cal timescale of the cavity mode dynamics, dictated by 
k = Db(a; a )|77(a; a )| 2 with D},(uj) the density of outside 
modes, is much longer than the inverse of the bandwidth 
of Db(uj)\r)(uj)\ 2 . Typically, for a power-law density of 
states, we may take this bandwidth to be oj a , so in fact 
we assume here k -C 0J a . Then, a can be taken out of the 
integral and by averaging, assuming (bj(0)) = 0, we get, 



ig(fi). 



(C3) 



The typical timescale for changes in the atomic operators 
is dictated by the decay rate to the perpendicular modes, 
which is similar to that in free-space, 70, and by g. By 
assuming the bad-cavity limit, i.e. k 3> ff,7o, and noting 
that (ft) < O(l), we have (a) ~ — f (a), and thus obtain 

(o)(t)«<o)(0)e-3 t . (C4) 
Finally, by the quantum regression theorem [19] we have 

. (d(t)tf(t + T )) = (a(f)at(t))e-* w = K + i) e -tH 

(at(t)a(t + T)) = (a*{t)a{t))e-^ T \ =n a e-i\ T \ (C5) 



where in the last step we assumed stationarity of the 
correlations, and took the thermal state as the stationary 
state, with n a the Planck distribution for uj a . 



Appendix D: Dipole matrix element calculation 

We would like to calculate the dipole matrix element 
d between the states |e) = \n,n — l,n — 1) and \g) = 
\n,n — 2, n — 2), with n = 51. Let us first derive an 
expression for a general n. Clearly the z component of 
the dipole vanishes since Am — 1, so let us calculate 
d x = er sin 9 cos <j), with the usual convention of spheri- 
cal coordinates. In terms of the Hydrogen wavefunctions 
Rni(r)Yi m (£l), where O is the solid angle, the matrix el- 
ement is written, 



d 

e 

R 



RA 



^ 3 K,„-i(^)^«,„- 2 (r) 



A = / dm. 



n— l,n— 1 



(fi) sin 6 cos 0y„_ 2 ,„_ 2 (fi)(Dl) 
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radial functions [24], we find 



Rn,n-l(r) 



1 



first noting that sin 9 cos 4> — J ^(Y^-i — and then 



We begin with the radial part R. Recalling the Hydrogen where T(x) is the Gamma function. For n = 51 we get 

R = — 768.815a . The angular part A is calculated by 

first ni 
using 

dftYi m Yi iq Yi' m ' = 



\nao J 2n(2n — 1)! 

ri-1 



e nao 



2r Y 
na / 



+ (2n-2)(^ 



n— 2 



Vna o y 2n(2n-2)! 



e nao 



nao / 
/ 2r 



3(2^ + 1) r 

47r(2Z + 1) ! ' 0;i '' 1 ' ' ,\,m' ,qi 

(D4) 

where C 'j,M-j 1 ,j 2 .m 1 ,m 2 are the Clebsch-Gordan coeffi- 
cients. We thus obtain, 



(D2) 



Performing the integrations in R, we obtain 



R 



A = 



2n- 3 



Cfl-l,0;li-2,l,0,0 



(2n- l)!(2n- 2)! 2 



(n - l)r(2n + 1) - 2 F ( 2n + 2 ) 



2(2n- 1) 

X [Cn-l,ri-l;n-2,l.n-2,-l _ Cra-l,n-l;n-2,l,n-2,l] • 

(D5) 

For n — 51 we get A = 0.497, so finally we find d = 
oq, (D3) -382.101ea . 



[1] P. D. Drummond and Z. Ficek (Eds.), Quantum Squeez- 
ing (Springer- Verlag, Berlin Heidelberg, 2004). 

[2] B. J. Dalton, Z. Ficek and S. Swain, J. Mod. Opt. 46, 
379 (1999). 

[3] C. W. Gardiner, Phys. Rev. Lett. 56, 1917 (1986). 

[4] H. J. Carmichael, A. S. Lane and D. F. Walls, Phys. Rev. 

Lett. 58, 2539 (1987). 
[5] H. T. Dung and N. Q. Khanh, J. Mod. Opt. 44, 1497 

(1997). 

[6] E. Shahmoon, S. Levit and R. Ozeri, Phys. Rev. A 80, 
033803 (2009). 

[7] A. S. Parkins, P. Zoller, and H. J. Carmichael, Phys. Rev. 

A 48, 758 (1993). 
[8] N. Liitkenhaus, J. I. Cirac and P. Zoller, Phys. Rev. A 

57, 548 (1998). 
[9] C. A. Muschik, E. S. Polzik, and J. I. Cirac, Phys. Rev. 

A 83, 052312 (2011). 
[10] G. Nikoghosyan and M. Fleischhauer, Phys. Rev. Lett. 

103, 163603 (2009). 
[11] S. G. Clark and A. S. Parkins, Phys. Rev. Lett. 90, 

047905 (2003). 

[12] R. Tanas, J. Opt. B: Quantum Semiclass. Opt. 4, S142- 
S152 (2002). 

[13] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-Photon Interactions: Basic Processes and Applica- 
tions, ( WILEY- VCH, 2004). 

[14] S. Maniscalco, J. Phlo and K.-A. Suominen, Eur. Phys. 
J. D 55, 181 (2009). 

[15] A. G. Kofman and G. Kurizki, Phys. Rev. A 54, R3750 
(1996); A. G. Kofman and G. Kurizki, Phys. Rev. Lett. 
87, 270405 (2001); A. G. Kofman and G. Kurizki, Phys. 



Rev. Lett. 93, 130406 (2004); A. G. Kofman and G. Kur- 
izki, Nature 405, 546 (2000). 

[16] G. Gordon, N. Erez and G. Kurizki, J. Phys. B 40, 
S75(2007); J. Clausen, G. Bensky, and G. Kurizki, Phys. 
Rev. Lett. 104, 040401 (2010); D. Vitali and P. Tombesi, 
Phys. Rev. A 59, 4178 (1999); D. Vitali and P. Tombesi, 
Phys. Rev. A 65, 012305 (2001). 

[17] N. Erez, G. Gordon, M. Nest and G. Kurizki, Nature 
452, 724-727 (2008); M. Kolaf, D. Gelbwaser-Klimovsky, 
R. Alicki and G. Kurizki, Phys. Rev. Lett. 109, 090601 
(2012). 

[18] D. F. Walls and G. J. Milburn, Quantum Optics, 
(Springer, 1995). 

[19] H. J. Carmichael, Statistical Methods in Quantum Optics 
1 , (Springer, 1998). 

[20] B. R. Mollow, Phys. Rev. A 12, 1919 (1975). 

[21] M.O. Scully and M. S. Zubairy, Quantum Optics (Cam- 
bridge University Press, Cambridge, England, 1997). 

[22] G. Gordon and G. Kurizki, Phys. Rev. Lett. 97, 110503 
(2006); Y. Khodorkovsky, G. Kurizki, and A. Vardi, 
Phys. Rev. Lett. 100, 220403 (2008); Y. Khodor- 
kovsky, G. Kurizki, and A. Vardi, Phys. Rev. A 80, 
023609 (2009); Nir Bar-Gill, Gershon Kurizki, Markus 
Oberthaler, and Nir Davidson, Phys. Rev. A 80, 053613 
(2009); I.E. Mazets and G. Kurizki, J. Phys. B 40, 
F105(2007). 

[23] A. Blais, R. S. Huang, A. Wallraff, S. M. Girvin and R. 

J. Schoelkopf, hys. Rev. A 69, 062320 (2004). 
[24] D. J. Griffiths, Introduction to Quantum Mechanics, 

(Prentice Hall, 1995). 



